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Abstract. In the search for appropriate discretizations of surface theory it is 
crucial to preserve such fundamental properties of surfaces as their invari- 
ance with respect to transformation groups. We discuss discretizations based 
on Mobius invariant building blocks such as circles and spheres. Concrete 
problems considered in these lectures include the Willmore energy as well 
as conformal and curvature line parametrizations of surfaces. In particular 
we discuss geometric properties of a recently found discrete Willmore en- 
ergy. The convergence to the smooth Willmore functional is shown for special 
refinements of triangulations originating from a curvature line parametriza- 
tion of a surface. Further we treat special classes of discrete surfaces such as 
isothermic and minimal. The construction of these surfaces is based on the 
theory of circle patterns, in particular on their variational description. 



1. Why from Circles? 

The theory of polyhedral surfaces aims to develop discrete equivalents of the geo- 
metric notions and methods of smooth surface theory. The latter appears then as 
a limit of refinements of the discretization. Current interest in this field derives 
not only from its importance in pure mathematics but also from its relevance for 
other fields like computer graphics. 

One may suggest many different reasonable discretizations with the same 
smooth limit. Which one is the best? In the search for appropriate discretiza- 
tions, it is crucial to preserve the fundamental properties of surfaces. A natural 
mathematical discretization principle is the invariance with respect to transfor- 
mation groups. A trivial example of this principle is the invariance of the the- 
ory with respect to Euclidean motions. A less trivial but well-known example is 
the discrete analog for the local Gaussian curvature defined as the angle defect 
G{v) = 27r — ^a^, at a vertex d of a polyhedral surface. Here the ai are the 
angles of all polygonal faces (see Fig. Isl) of the surface at vertex v. The discrete 
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Figure 1. Discrete surfaces made from circles: general simplicial 
surface and a discrete minimal Enneper surface. 

Gaussian curvature G{v) defined in this way is preserved under isometrics, which 
is a discrete version of the theoremum egregium of Gauss. 

In these lectures, we focus on surface geometries invariant under Mobius 
transformations. Recall that Mobius transformations form a finite-dimensional Lie 
group generated by inversions in spheres (see Fig. [2|. Mobius transformations 




Figure 2. Inversion B i-^ C in a sphere, | AB \\ AC |= i?^. A 
sphere and a torus of revolution and their inversions in a sphere: 
spheres are mapped to spheres. 

can be also thought as compositions of translations, rotations, homothcties and 
inversions in spheres. Alternatively, in dimensions n > 3 Mobius transformations 
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can be characterized as conformal transformations: Due to Liouville's theorem any 
conformal mapping F : U ^ V between two open subsets U,V C M" , n > 3 is a 
Mobius transformation. 

Many important geometric notions and properties are known to be preserved 
by Mobius transformations. The hst includes in particular: 

• spheres of any dimension, in particular circles (planes and straight lines are 
treated as infinite spheres and circles), 

• intersection angles between spheres (and circles), 

• curvature line parametrization, 

• conformal parametrization, 

• isothermic parametrization ( conformal curvature line parametrization), 

• the Willmore functional (see Section [2| . 

For discretization of Mobius-invariant notions it is natural to use Mobius- 
invariant building blocks. This observation leads us to the conclusion that the 
discrete conformal or curvature line parametrizations of surfaces and the discrete 
Willmore functional should be formulated in terms of circles and spheres. 

2. Discrete Willmore Energy 

The Willmore functional [39j for a smooth surface S in 3-dimensional Euclidean 
space is 

W{S) = J [ {ki- k2fdA = [ H'^dA - [ KdA. 

Here dA is the area element, fci and ^2 the principal curvatures, H = ^(fci + ^2) 
the mean curvature, and K = kik2 the Gaussian curvature of the surface. 
Let us mention two important properties of the Willmore energy: 

• yV(S') > and W{S) = if and only if 5 is a round sphere. 

• W{S) (and the integrand {ki — k2)'^dA) is Mobius invariant [11 [55] 

Whereas the first claim almost immediately follows from the definition, the second 
is a non-trivial property. Observe that for closed surfaces W{S) and H^dA differ 
by a topological invariant / KdA = 27rx(>5'). We prefer the definition of yV(5) with 
a Mobius invariant integrand. 

Observe that minimization of the Willmore energy W seeks to make the 
surface "as round as possible" . This property and the Mobius invariance are two 
principal points of the geometric discretization of the Willmore energy suggested 
in In this section we present the main results of [3j with complete derivations, 
some of which were omitted there. 

2.1. Discrete Willmore functional for simplicial surfaces 

Let 5* be a simplicial surface in 3-dimensional Euclidean space with vertex set V , 
edges E and (triangular) faces F. We define the discrete Willmore energy of S 
using the circumcircles of its faces. Each (internal) edge e € E is incident to two 
triangles. A consistent orientation of the triangles naturally induces an orientation 
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of the corresponding circumcircles. Let /3(e) be the external intersection angle of 
the circumcircles of the triangles sharing e, meaning the angle between the tangent 
vectors of the oriented circumcircles (at either intersection point). 

Definition 2.1. The local discrete Willmore energy at a vertex v is the sum 

W{v) = ^/3(e) - 27r. 

over all edges incident to v. The discrete Willmore energy of a compact simplicial 
surface S without boundary is the sum over all vertices 

m^) = ^E^(^) = E/3(e)-^i^i- 

Here \V\ is the number of vertices of S. 




Figure 3. Definition of discrete Willmore energy 

Figure |3]presents two neighboring circles with their external intersection angle 
f3i as well as a view "from the top" at a vertex v showing all n circumcircles passing 
through V with the corresponding intersection angles /3i, . . . , /3„. For simplicity we 
will consider only simplicial surfaces without boundary. 

The energy W{S) is obviously invariant with respect to Mobius transforma- 
tions. 

The star S{v) of the vertex v is the subcomplex of S consisting of the triangles 
incident with v. The vertices of S(y) are v and all its neighbors. We call S{v) convex 
if for each of its faces / S F{S{v)) the star S{v) lies to one side of the plane of / 
and strictly convex if the intersection of S (v) with the plane of / is / itself. 

Proposition 2.2. The conformal energy is nonnegative 

W{v) > 0, 

and vanishes if and only if the star S{v) is convex and all its vertices lie on a 
common sphere. 

The proof of this proposition is based on an elementary lemma. 
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Lemma 2.3. Let V he a (not necessarily planar) n-gon with external angles Pi. 
Choose a point P and connect it to all vertices of V . Let ai be the angles of the 
triangles at the tip P of the pyramid thus obtained (see Figure^. Then 

n n 

'i=l 1=1 

and equality holds if and only ifV is planar and convex and the vertex P lies inside 




Figure 4. Proof of Lemma 12.31 

Proof. Denote by 7^ and 5i the angles of the triangles at the vertices oi V, as in 
Figure |4] The claim of Lemma |2.3| follows from adding over all i ~ 1 , . . . , rt the 
two obvious relations 

(ii+i > TT - (7j+i + Si) 
= 7r-(7j + (5i). 

All inequalities become equalities only in the case when V is planar, convex and 
contains P. □ 

For P in the convex hull of P we have Q^i ^ 27r. As a corollary we obtain 
a polygonal version of Fenchel's theorem ^TE\ : 

Corollary 2.4. 

n 

Proof of Proposition \2.S\ The claim of Proposition |2.2| is invariant with respect 
to Mobius transformations. Applying a Mobius transformation M that maps the 
vertex v to infinity, M{v) — 00, we make all circles passing through v into straight 
lines and arrive at the geometry shown in Figure |4] with P = M{oo). Now the 
result follows immediately from Corollary |2.4[ 

Theorem 2.5. Let S be a compact simplicial surface without boundary. Then 

W{S) > 0, 

and equality holds if and only if S is a convex polyhedron inscribed in a sphere, 
i.e. a Delaunay triangulation of a sphere. 
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Proof. Only the second statement needs to be proven. By Proposition |2.2| the 
equality W{S) = implies that the star of each vertex of S is convex (but not 
necessarily strictly convex). Deleting the edges that separate triangles lying in a 
common plane, one obtains a polyhedral surface Sp with circular faces and all 
strictly convex vertices and edges. Proposition |2.2| implies that for every vertex v 
there exists a sphere with all vertices of the star S{v) lying on it. For any edge 
(i'i,W2) of Sp two neighboring spheres Sy-^ and Sy^ share two different circles of 
their common faces. This implies Sy-^ — Sy^ and finally the coincidence of all the 
spheres Sy. □ 

2.2. Noninscribable polyhedra 

The minimization of the conformal energy for simplicial spheres is related to a 
classical result of Steinitz [37], who showed that there exist abstract simplicial 
3-polytopes without geometric realizations as convex polytopes with all vertices 
on a common sphere. We call these combinatorial types noninscribable. 

Let 5* be a simplicial sphere with vertices colored in black and white. Denote 
the sets of white and black vertices by Vyj and 14 respectively, V = Vyj^Vb- Assume 
that there are no edges connecting two white vertices and denote the sets of the 
edges connecting white and black vertices and two black vertices by Ey,b and Ebb 
respectively, E = Ey,b U Ebb- The sum of the local discrete Willmore energies over 
all white vertices can be represented as 



^ W{v) = J2 - 27r|K, 



Its nonnegativity yields J^eeE ^ f^^^^ — For the discrete Willmore energy 

of S this implies 

Wis) = /3(e) + E ^(^) - ^(l^-l + 1^"!) ^ '^(l^-l - l^"!)- (2.1) 

The equality here holds if and only if /3{e) — for all e G Ebb and the star of any 
white vertices is convex, with vertices lying on a common sphere. We come to the 
conclusion that the polyhedra of this combinatorial type with > |Vb| have 
positive Willmore energy and thus cannot be realized as convex polyhedra all of 
whose vertices belong to a sphere. These are exactly the noninscribable examples 
of Steinitz pi] . 

One such example is presented in Figure [5] Here the centers of the edges of 
the tetrahedron are black and all other vertices are white, so \Vyj\ — 8, \ Vb\ — 6. 



The estimate (2.11 implies that the discrete Willmore energy of any polyhedron 
of this type is at least 27r. The polyhedra with energy equal to 27r are constructed 
as follows. Take a tetrahedron, color its vertices white and chose one black vertex 
per edge. Draw circles through each white vertex and its two black neighbors. We 



get three circles on each face. Due to Miquel's theorem (see Fig. 10) these three 



circles intersect at one point. Color this new vertex white. Connect it by edges 
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to all black vertices of the triangle and connect pairwise the black vertices of the 
original faces of the tetrahedron. The constructed polyhedron has W — 2tt. 





Figure 5. Discrete Willmore spheres of inscribable {W — 0) and 
non-inscribable {W > 0) types. 

To construct further polyhedra with > |Vf,|, take a polyhedron P whose 
number of faces is greater then the number of vertices |F| > \V\. Color all the 
vertices black, add white vertices at the faces and connect them to all black vertices 
of a face. This yields a polyhedron with \Vw\ = \F\ > \Vb\ — \V\. Hodgson, Rivin 
and Smith [24J have found a characterization of inscribable combinatorial types, 
based on a transfer to the Klein model of hyperbolic 3-space. Their method is 
related to the methods of construction of discrete minimal surfaces in Section [S] 

The example in Figure [5] (right) is one of the few for which the minimum 
of the discrete Willmore energy can be found by elementary methods. Generally 
this is a very appealing (but probably difficult) problem of discrete differential 
geometry (see the discussion in [3J). 

Complete understanding of noninscribable simplicial spheres is an interest- 
ing mathematical problem. However the phenomenon of existence of such spheres 
might be seen as a problem in using of the discrete Willmore functional for ap- 
plications in computer graphics, such as the fairing of surfaces. Fortunately the 
problem disappears after just one refinement step: all simplicial spheres become 
inscribable. Let S be an abstract simplicial sphere. Define its refinement Sr as 
follows: split every edge of S in two by inserting additional vertices, and connect 
these new vertices sharing a face of S by additional edges (1^4 refinement, see 
Figure |7] (left)). 

Proposition 2.6. The refined simplicial sphere Sr is inscribable, and thus there 
exists a polyhedron Sr with the combinatorics o/ Sr and W{Sii) = 0. 



Proof. Koebe's theorem (see Theorem 5.3 Section p|) states that every abstract 



simplicial sphere S can be realized as a convex polyhedron S all of whose edges 
touch a common sphere . Starting with this realization S it is easy to construct 
a geometric realization Sr of the refinement Sr inscribed in S^. Indeed, choose 
the touching points of the edges of S with S'^ as the additional vertices of Sr and 
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project the original vertices of S (which lie outside of the sphere 5^) to S^. One 
obtains a convex simplicial polyhedron Sji inscribed in 5*^. □ 

2.3. Computation of the energy 

For derivation of some formulas it will be convenient to use the language of quater- 
nions. Let {l,i,j,k} be the standard basis 

ij = k, ij = k, ij = k, ii = jj = kk = 1 

of the quaternion algebra H. A quaternion q = go 1 + 91 i + '?2j + 13^ is decomposed 
in its real part Re q := qq € M. and imaginary part Im q := qii + q2i + f73k €E Im H. 
The absolute value of q is \q\ := f/o + 9i + 92 + 93- 

We identify vectors in M"^ with imaginary quaternions 

V = {vi,V2, W3) G < > V = vi'i + V2j + v-^k e Im IHI 

and do not distinguish them in our notation. For the quaternionic product this 
implies 

vw ~ —{v,w) + V X w, (2-2) 
where {v, w) and v x w are the scalar and vector products in M^. 

Definition 2.7. Let xi,X2,xz,X4 G M.^ = Im H be points in 3-dimensional Euclidean 
space. The quaternion 

q{xi,X2,X3,X4) := {xi - X2){X2 - X3)~^{X3 - X4){X4 - Xi)~^ 

is called the cross-ratio of xi,X2,X3,X4. 

The cross-ratio is quite useful due to its Mobius properties: 

Lemma 2.8. The absolute value and real part of the cross-ratio q{xi, X2, x^^ X4) (or 
equivalently \q\ and |Im q\) are preserved by Mobius transformations. The quadri- 
lateral xi,X2,X3,X4 is circular if and only if q(xi, X2, X3, X4) € M. 

Consider two triangles with a common edge. Let a, b,c,d E M.^ be their other 
edges, oriented as in Figj6] 




Figure 6. Formula for the angle between circumcircles 
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Proposition 2.9. The external angle [3 G [0,7r] between the circumcircles of the 
triangles in Figure\S^is given by any of the equivalent formulas: 

Kc q Re (abed) 

= -Ji\=-^^b^ = 
^ (q, c) {b, d) - (g, b) (c, d) ~ {b, c) {d, a) 
\a\\b\\c\\d\ ■ 

Here q = ab^^cd^^ is the cross-ratio of the quadrilateral. 

Proof. Since Re q, \q\ and f3 are Mobius- invariant, it is enough to prove the first 
formula for the planar case a,b,c,d G C, mapping all four vertices to a plane 
by a Mobius transformation. In this case q becomes the classical complex cross- 
ratio. Considering the arguments a,b,c^d € C one easily arrives at /3 = tt — arg q. 
The second representation follows from the identity b^^ = —b/\b\ for imaginary 
quaternions. Finally applying (|2.2| we obtain 



Re (abed) = (a, b) (c, d) — {a x b,c x d) = 
(a, b) {c,d) + (b, c) {d, a) ~ (a, c) {b,d). 

□ 



2.4. Smooth limit 

The discrete energy W is not only a discrete analogue of the Willmore energy. In 
this section we show that it approximates the smooth Willmore energy, although 
the smooth limit is very sensitive to the refinement method and should be chosen 
in a special way. We consider a special infinitesimal triangulation which can be 
obtained in the limit of 1 ^ 4 refinements (see Fig. [t] (left)) of a triangulation 
of a smooth surface. Intuitively it is clear that in the limit one has a regular 
triangulation such that almost every vertex is of valence 6 and neighboring triangles 
are congruent up to sufficiently high order of e (e is the order of the distances 
between neighboring vertices). 




Figure 7. Smooth limit of the discrete Willmore energy. Left: 
The 1 — > 4 refinement. Middle: An infinitesimal hexagon in the 
parameter plane with a (horizontal) curvature line. Right: The 
/3-angle corresponding to two neighboring triangles in M.^. 
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We start with a comparison of the discrete and smooth WiUmore energies for 
an important modeUing example. Consider a neighborhood of a vertex v Cz S, and 
represent the smooth surface locally as a graph over the tangent plane at v: 



9 {x,y)^ f{x,y)^ (x,y,^{hx^ + k2y^) + oix^ + y^)^ eR^, {x,y) 



0. 



Here x,y are the curvature directions and ki, k2 are the principal curvatures at v. 
Let the vertices (0, 0), a = (ai, 02) and b = (61, 62) in the parameter plane form an 
acute triangle. Consider the infinitesimal hexagon with vertices ea, eb, ec, —ea, —eb, — ec, 
(see Figure I?] (middle)), with b = a+c. The coordinates of the corresponding points 
on the smooth surface are 

/(±ea) = e(±ai,±a2,era + o(e)), 
/(±ec) = e(±Ci,±C2,erc + o(e)), 

f{±eb) = {f{±ea) + f{±ec))+e^R, R ^ {0,0,r + o{e)), 

where 

ra^^{kial + k2al), ^ ^{kicj + k2cl), r = {kiUiCi + k2a2C2) 

and a = (ai, 02), c = (ci, C2). 

We will compare the discrete Willmore energy W of the simplicial surface 
comprised by the vertices /(ea), . . . , f{—ec) of the hexagonal star with the classical 
Willmore energy W of the corresponding part of the smooth surface S. Some 
computations are required for this. Denote by eA = f{ea), eB = f{eb), eC = f{ec) 
the vertices of two corresponding triangles (as in Figure [t] (right)), and also by \a\ 
the length of a and by (a, c) = oiCi + 0202 the corresponding scalar product. 

Lemma 2.10. The external angle /3(e) between the circumcircles of the triangles 
with the vertices {0,A,B) and {0,B,C) (as in Figure^ (right)) is given by 



Here /3(0) is the external angle of the circumcircles of the triangles (0, a, b) and 
(0, 6, c) in the plane, and 

g = |apr,(r + r,) + |cpr,(r + ra) + y(|ap + |cn, 

h = |a|Ve(r + rc) + |c|Va(r + ra) - (a,c)(r + 2ra)(r + 2rc). 



Proof Formula ( 2.3 1 with a = -C,b = A,c = C + eR,d = -A - eR yields 



{C\C + eR){A,A + eR) - {A,C){A + eR,C + eR) - {A,C + eR){A + eR,C) 
""^"^ ~ \A\\C\\A + eR\\C + eR\ ' 

where \A\ is the length of A. Substituting the expressions for A,C,R we see that 
the term of order e of the numerator vanishes, and we obtain for the numerator 

\a\^\c\^-2{a,c)^ + e^h + o{e^). 
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For the terms in the denominator we get 

1^1 = l«l (l + + 0(6^)) , \A + eR\ = \a\ [l + ^^^^^^^ + o{e^) 

and similar expression for \C\ and \C + eR\. Substituting this to the formula for 
cos /3 we obtain 

Observe that this formula can be read as 

cos /3(e) = cos/3(0) + . |„| .„ (/i - gcos/3(0)) + o(e^ 




which implies the asymptotic (2.4 1. □ 

The term w{h) is in fact the part of the discrete Willmore energy of the vertex 
V coming from the edge b. Indeed the sum of the angles /3(0) over all 6 edges 
meeting at v is 2n. Denote by w(a) and w{c) the parts of the discrete Willmore 
energy corresponding to the edges a and c. Observe that for the opposite edges 
(for example a and —a) the terms w coincide. Denote We{v) the discrete Willmore 
energy of the simplicial hexagon we consider. We have 

W^{v) = (w(a) + w{h) + w{c)) + o(e^). 

On the other hand the part of the classical Willmore functional corresponding to 
the vertex v is 

W,(z;)= J(fci-fc2)25 + 0(e2), 

where the area S is one third of the area of the hexagon or, equivalently, twice the 
area of one of the triangles in the parameter domain 

S — e^\a\\c\ sm"f. 

Here 7 is the angle between the vectors a and c. An elementary geometric consid- 
eration implies 

/3(0) = 27~7r. (2.5) 

We are interested in the quotient W^/Wc which is obviously scale invariant. Let us 
normalize |a| = 1 and parametrize the triangles by the angles between the edges 
and by the angle to the curvature line (see Fig. |7] (middle)). 

(01,02) = (cos 01, sin 01 ), (2.6) 

, . , sin 02 / , , , s sin02 . / , , , 

(Cl,C2) = {- — COS(0i + 02 + 03), ^ — Sm(0i + 02 + 03)). 

sm 03 sm 03 

The moduli space of the regular lattices of acute triangles is described as follows 

$= {0= (01,02,03) e M^|0 < 01 < |, < 02 < |, < 03 < |, I < 02+03}- 
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Proposition 2.11. The limit of the quotient of the discrete and smooth Willmore 
energies 

is independent of the curvatures of the surface and depends on the geometry of the 
triangulation only. It is 

0(6) - 1 - ('^Q'^^'^i cos 03 + cos(2(/)i + 2(j)2 + (j)3))^ + (sin 201 cos^s)^ 

4 cos (j>2 cos 03 COS(02 + 03) 

and we have Q > I. The infimum inf$ Q(0) = 1 corresponds to one the cases when 
two of the three lattice vectors a,b,c are in the principal curvature directions 

• 01 = 0,02 + 03 -> f, 

• 01=0,02^1, 

• 01+02= 1,03 -> f. 

Proof, is based on a direct but rather involved computation. We have used the 
Mathematica computer algebra system for some of the computations. Introduce 

4w 



■ (fci-fc2)'^' 

This gives in particular 

^ fe + g(2cos^7-l) ^ h + 5(2^°ifJ - 1) 



(fci -fc2)2|a|3|c|3cos7sin2 7 (fci - fc2)2(a, c)(|a|2|c|2 - (a, c)2) ' 

Here we have used the relation (2.51 between /3{0) and 7. In the sum over the 
edges Q = w{a) + w^b) + u;(c) the curvatures ki,k2 disappear and we get Q in 
terms of the coordinates of a and c: 

Q = 2{{alcl + alcl){aici + a2C2) + alcl{al + cl) + alcl{al + cf) + 
2aia2CiC2((ai + Ci)^ + (02 + ca)^))/ 

((aici + a2C2)(ai(ai + ci) + 02(02 + C2))((ai + ci)ci + (02 + C2)c2)) ■ 
Substituting the angle representation ( |2.6[ ) we obtain 

^ _ sin 201 sin 2(0i + 02) + 2 cos 02 sin(20i + 02) sin 2(0i + 02 + 03) 

4 cos 02 cos 03 cos ( 02 + 03 ) 

One can check that this formula is equivalent to (2.7 1. Since the denominator in 
(2.7) on the space $ is always negative we have Q > 1. The identity Q = 1 holds 
only if both terms in the nominator of (2.7 1 vanish. This leads exactly to the 
cases indicated in the proposition when the lattice vectors are directed along the 
curvature lines. Indeed the vanishing of the second term in the nominator implies 
either 0i = or 03 ^ ^ . Vanishing of the first term in the nominator with 0i = 
implies 02 ^ f or 02 + 03 — > ^ . Similarly in the limit 03 — > § the vanishing of 

(cos 201 cos 03 + cos(20i + 202 + 03))^/cOS03 

implies 0i + 02 = f • One can check that in all these cases (3(0) — > 1. □ 
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Note that for the infinitesimal equilateral triangular lattice <j)2 ~ <f>3 — ^ the 
result is independent of the orientation (f)i with respect to the curvature directions, 
and the discrete Willmore energy is in the limit Q=3/2 times larger than the 
smooth one. 

Finally, we come to the following conclusion. 

Theorem 2.12. Let S be a smooth surface with Willmore energy yV(5) . Consider 
a simplicial surface such that its vertices lie on S and are of degree 6, the 
distances between the neighboring vertices are of order e, and the neighboring tri- 
angles of meeting at a vertex are congruent up to order (i.e. the lengths of 
the corresponding edges differ by terms of order at most e*), and they build ele- 
mentary hexagons the lengths of whose opposite edges differ by terms of order at 
most e^. Then the limit of the discrete Willmore energy is bounded from below by 
the classical Willmore energy 

limW(S',) > >V(5). (2.8) 



Moreover equality in (2.8) holds if is a regular triangulation of an infinitesimal 
curvature line net of S , i.e. the vertices of are at the vertices of a curvature 
line net of S. 

Proof. Consider an elementary hexagon of S^. Its projection to the tangent plane 
of the central vertex is a hexagon which can be obtained from the modelling 



one considered in Proposition 2.11 by a perturbation of vertices of order o{e^). 
Such perturbations contribute to the terms of order o(e^) of the discrete Willmore 
energy. The latter are irrelevant for the considerations of Proposition |2.11[ □ 

Possibly minimization of the discrete Willmore energy with the vertices con- 
strained to lie on S could be used for computation of a curvature- line net. 

2.5. Bending energy for simplicial surfaces 

An accurate model for bending of discrete surfaces is important for modeling in 
computer graphics. 

The bending energy of smooth thin shells is given by the integral |19j 



where Hq and H are the mean curvatures of the original and deformed surface 
respectively. For ifg = it reduces to the Willmore energy. 

To derive the bending energy for simplicial surfaces let us consider the limit of 
fine triangulations, where the angles between the normals of neighboring triangles 
become small. Consider an isometric deformation of two adjacent triangles. Let 9 
be the external dihedral angle of the edge e, or, equivalently, the angle between 
the normals of these triangles (see Figure [s]) and f3{9) the external intersection 
angle between the circumcircles of the triangles (see Figure [3]) as a function of 8. 
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Proposition 2.13. Assume that the circumcenters of two adjacent triangles do not 
coincide. Then in the limit of small angles 9 —^ the angle (3 between the circles 
behaves as follows: 

Here I is the length of the edge and L =/= is the distance between the centers of 
the circles. 

Proof. Let us introduce the orthogonal coordinate system with the origin at the 
middle point of the common edge e, the first basis vector directed along e, and 
the third basis vector orthogonal to the left triangle. Denote by Xi,X2 the centers 
of the circumcircles of the triangles and by , X4 the end points of the com- 
mon edge (see Fig. |8|. The coordinates of these points are Xi = (0, —li,0),X2 = 
(0,/2COs6',Z2sin6'),X3 = {h,0,0),,X4 = (-?3,0,0). Here 2I3 is the length of the 
edge e, and li and I2 are the distances from its middle point to the centers of the cir- 
cumcirlces (for acute triangles). The unit normals to the triangles are Ni = (0, 0, 1) 
and N2 — (0, — sin 61, cos 0). The angle /3 between the circumcircles intersecting at 
the point X4 is equal to the angle between the vectors A = Ni x {X4 — Xi) and 
B = N2 X {X4 — X2). The coordinates of these vectors are A = (—^1,-/3,0), 
B = {I2, —I3 cos 0, —I3 sin 6). This implies for the angle 

^/oN l\ cos 9 — UU 

cos/3(6') = ^ (2.9) 

rir2 

where = \Ar+^! i = 1, 2 are the radii of the corresponding circumcircles. Thus 
I3{9) is an even function, in particular (3{9) = /3(0) + B9^ + o{9^). Differentiating 



(2.9) by 9^ we obtain 



B = 



Also formula (2.9) yields 



2rir'2 sin /3(0) 



sin/3(0) 

rir2 
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where L — |Zi + ^2! is the distance between the centers of the circles. Finally 
combining these formulas we obtain B = 1^/{2L). □ 

This proposition motivates us to define the bending energy of simplicial sur- 
faces as 

eeE 

For discrete thin-shells this bending energy was suggested and analyzed by Grin- 
spun et al. [201 Us] ■ The distance between the barycenters was used for L in the 
energy expression, and possible advantages in using circumcenters were indicated. 
Numerical experiments demonstrate good qualitative simulation of real processes. 

Further applications of the discrete Willmore energy in particular for surface 
restoration, geometry denoising, and smooth filling of a hole can be found in [8j . 

3. Circular Nets as Discrete Curvature Lines 

Simplicial surfaces as studied in the previous section are too unstructured for 
analytical investigation. An important tool in the theory of smooth surfaces is 
the introduction of (special) parametrizations of a surface. Natural analogues of 
parametrized surfaces are quadrilateral surfaces, i.e. discrete surfaces made from 
(not necessarily planar) quadrilaterals. The strips of quadrilaterals obtained by 
gluing quadrilaterals along opposite edges can be considered as coordinate lines 
on the quadrilateral surface. 

We start with a combinatorial description of the discrete surfaces under con- 
sideration. 

Definition 3.1. A cellular decomposition 2? of a two-dimensional manifold (with 
boundary) is called a quad-graph if the cells have four sides each. 

A quadrilateral surface is a mapping / of a quad-graph to . The mapping 
/ is given just by the values at the vertices of V, and vertices, edges and faces of 
the quad-graph and of the quadrilateral surface correspond. Quadrilateral surfaces 
with planar faces were suggested by Sauer |32| as discrete analogs of conjugate nets 
on smooth surfaces. The latter are the mappings {x,y) ^ f{x,y) € such that 
the mixed derivative fxy is tangent to the surface. 

Definition 3.2. A quadrilateral surface / : I? ^ all faces of which are circular 
(i.e. the four vertices of each face lie on a common circle) is called a circular net 
(or discrete orthogonal net). 

Circular nets as discrete analogues of curvature line parametirzed surfaces 
were mentioned by Martin, de Pont, Sharrock and Nutbourne [23 130] • The cur- 
vature lines on smooth surfaces continue through any point. Keeping in mind the 
analogy to the curvature line parametrized surfaces one may in addition require 
that all vertices of a circular net are of even degree. 
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A smooth conjugate net / : Z? — )■ M"^ is a curvature line parametrization if 
and only if it is orthogonal. The bisectors of the diagonals of a circular quadrilat- 
eral intersect orthogonally (see Figure |9]) and can be interpreted [TTI as discrete 
principal curvature directions. 




Figure 9. Principal curvature directions of a circular quadrilat- 
eral 



There are deep reasons to treat circular nets as a discrete curvature line 
parametrization. 

• The class of circular nets as well as the class of curvature line parametrized 
surfaces is invariant under Mobius transformations. 

• Take an infinitesimal quadrilateral (/(x, y), /(x-t-e), y), f{x+e), y+e), f{x, y+ 
e)) of a curvature line parametrized surface. A direct computation (see [11 J 
shows that in the limit e — > the imaginary part of its cross-ratio is of 
order e^. Note that circular quadrilaterals are characterized by having real 
cross-ratios. 

• For any smooth curvature line parametrized surface f : D M!^ there exists 
a family of discrete circular nets converging to /. Moreover the convergence 
is C°° , i.e. with all derivatives. The details can be found in |5]. 

One more argument in favor of Definition |3.2| is that circular nets satisfy the 
consistency principle, which has proven to be one of the organizing principles in 
discrete differential geometry |10j . The consistency principle singles out fundamen- 
tal geometries by the requirement that the geometry can be consistently extended 
to a combinatorial grid one dimension higher. The consistency of circular nets was 
shown by Doliwa and Santini I16j based on the following classical theorem. 

Theorem 3.3. (Miquel) Consider a combinatorial cube in M'^ with planar faces. 
Assume that three neighboring faces of the cube are circular. Then the remaining 
three faces are also circular. 

Equivalently, provided the four-tuples of black vertices corresponding to three 
neighboring faces of the cube lie on circles, the three circles determined by the 
triples of points corresponding to three remaining faces of the cube all intersect 



(at the white vertex in Figure 10 1. It is easy to see that all vertices of Miquel's 
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cube lie on a sphere. Mapping the vertex shared by the three original faces to 
infinity by a Mobius transformation, we obtain an equivalent planar version of 
Miquel's theorem. This version, also shown in Figure \T0\ can be proven by means 
of elementary geometry. 




Figure 10. Miquel's theorem: spherical and planar versions 



Finally note that circular nets are also treated as a discretization of triply or- 
thogonal coordinate systems. Triply orthogonal coordinate systems in K'^ are maps 
{x, y, z) 1-^ f{x, y, z) g E'^ from a subset of E.^ with mutually orthogonal fx, fy, fz- 
Due to the classical Dupin theorem, the level surfaces of a triply orthogonal coor- 
dinate system intersect along their common curvature lines. Accordingly, discrete 
triply orthogonal systems are defined as maps from (or a subset thereof) to M.^ 
with all elementary hexahedra lying on spheres Due to Miquel's theorem a dis- 
crete orthogonal system is uniquely determined by the circular nets corresponding 
to its coordinate two-planes (see [16] and [10]). 



4. Discrete Isothermic Surfaces 

In this section and in the following one, we investigate discrete analogs of special 
classes of surfaces obtained by imposing additional conditions in terms of circles 
and spheres. 

We start with minor combinatorial restrictions. Suppose that the vertices of 
V are colored black or white so that the two ends of each edge have different colors. 
Such a coloring is always possible for topological discs. To model the curvature 
lines, suppose also that the edges of a quad-graph V may consistently be labelled 
'-I-' and '— ', as in Figure [TT] (For this it is necessary that each vertex has an even 
number of edges). Let /o be a vertex of a circular net, fi, fs, ■ ■ ■ , f4N-i be its 



neighbors, and /2, fi,---, f^N its next-neighbors (see Figure 12 (left)). We call the 
vertex /o generic if it is not co-spherical with all its neighbors and a circular net 
f : V ^ M.^ generic if all its vertices are generic. 

Let / : I? ^ M"^ be a generic circular net such that every its vertex is co- 
spherical with all its next-neighbors. We will call the corresponding sphere central. 
For analytical description of this geometry let us map the vertex fo to infinity by 
a Mobius transformation A4{fo) — oo, and denote by Fi — A4{fi), i — 1, . . . ,4N 
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Figure 11. Labelling the edges of a discrete isothermic surface 



the Mobius images of fi. The points F2, F4, . . . , ^4^^ are obviously coplanar. The 
circles of the faces are mapped to straight lines. For the cross-ratios we get 

/ _p _r ( -c \ _ -^2A; — -F2/C + I Z2k + l 

Q(J0, J2k-1, J2k, J2k+l) — 7^ — , 

J'2k — ^2k~l Z2k~l 

where Z2k+i is the coordinate of F2k+i orthogonal to the plane V of F2, F4, . . . , F4n- 
(Note that since /o is generic none of the Zi vanishes). As a corollary we get for 
the product of all cross-ratios: 



Y\_ lifo, f2k-l, f2k, hk+l) — 1- 



(4.1) 



fc=i 



U h h 



1 


( 


1 


/o 



h h fs 




Figure 12. Central spheres of a discrete isothermic surface: com- 
binatorics (left), and the Mobius normalized picture for N = 2 
(right). 



Definition 4.1. A circular net / : I? — > M'^ satisfying condition (4.1 1 at each vertex 
is called a discrete isothermic surface. 

This definition was first suggested in ^ for the case of the combinatorial 
square grid 2? = Z^. In this case if the vertices are labelled by fm,n and the corre- 
sponding cross-ratios by q,n,„ := g(/m,«, /m+i,«, /m+i,n+i, /m,n+i), the condition 



(4.1) reads 

9m,n'7m+l,n+l '?m+l,n9m,n+l ■ 

Proposition 4.2. Every vertex fm,n of a discrete isothermic surface f : 1? ^ M? 
possesses a central sphere, i.e. the points fm,n, fm-i,n-i, frn+i,n-i, fm+i,n+i and 
fm-i,n+i are co-spherical. Moreover for generic circular maps f : Z"^ M.^ this 
property characterizes discrete isothermic surfaces. 
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Proof. Use the notation of Figure 12 with /o = fm.n, and the same argument 
with the Mobius transformation M which maps /o to oo. Consider the plane V 
determined by the points F2 , F4 and Fg . Let as above zj. be the coordinates of 



orthogonal to the plane V. Condition (4.1) yields 



Fs^Fr zr' 

This implies that the z-coordinate of the point Fg, vanishes, thus Fg E V. □ 

The property to be discrete isothermic is also SD-consistent, i.e. can be con- 
sistently imposed on all faces of a cube. This was proven first by Hertrich-Jeromin, 
Hoffmann and Pinkall [53] (see also [TU] for generalizations and modern treatment). 

An important subclass of discrete isothermic surfaces is given by the con- 
dition that all the faces are conformal squares, i.e. their cross ratio equals — 1. 
All conformal squares are Mobius equivalent, in particular equivalent to the stan- 
dard square. This is a direct discretization of the definition of smooth isothermic 
surfaces. The latter are immersions {x,y) 1-^ f{x,y) E M.^ satisfying 

||/x|| = WfyWJx -L fyjxy e span{/^,/J, (4.2) 

i.e. conformal curvature line parametrizations. Geometrically this definition means 
that the curvature lines divide the surface into infinitesimal squares. 



Figure 13. Right-angled kites are conformal squares 

The class of discrete isothermic surfaces is too general and the surfaces are 
not rigid enough. In particular one can show that the surface can vary preserving 
all its black vertices. In this case, one white vertex can be chosen arbitrarily [7]. 
Thus, we introduce a more rigid subclass. To motivate its definition, let us look at 
the problem of discretizing the class of conformal maps / : D ^ C for £) C C = R^. 
Conformal maps are characterized by the conditions 

\M^\fy\, f.±fy. (4.3) 

To define discrete conformal maps / : D ^ C, it is natural to impose these 
two conditions on two different sub-lattices (white and black) of Z^, i.e. to require 
that the edges meeting at a white vertex have equal length and the edges at a 
black vertex meet orthogonally. This discretization leads to the circle patterns 
with the combinatorics of the square grid introduced by Schramm [34] . Each circle 
intersects four neighboring circles orthogonally and the neighboring circles touch 



cyclically (Fig 14 (left)) 
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Figure 14. To the definition of discrete S-isothermic surfaces: 
orthogonal circle patterns as discrete conformal maps (left) and 
combinatorics of S-quad-graphs (right). 



The same properties imposed for quadrilateral surfaces with the combina- 
torics of the square grid f : 7? D "D ^ M.^ lead to an important subclass of discrete 
isothermic surfaces. Let us require for a discrete quadrilateral surface that: 

• the faces are orthogonal kites, 

• the edges meet at black vertices orthogonally (black vertices are at orthogonal 
corners of the kites), 

• the kites which do not share a common vertex are not coplanar (locality 
condition) . 

Observe that the orthogonality condition (at black vertices) implies that one 
pair of opposite edges meeting at a black vertex lies on a straight line. Together 
with the locality condition this implies that there are two kinds of white vertices, 
which we denote by © and (s). Each kite has white vertices of both types and the 
kites sharing a white vertex of the first kind © axe coplanar. 

These conditions imposed on the quad-graphs lead to S-quad-graphs and S- 
isothermic surfaces (the latter were introduced in [7] for the combinatorics of the 
square grid). 

Definition 4.3. An S-quad-graph I? is a quad-graph with black and two kinds of 
white vertices such that the two ends of each edge have different colors and each 



quadrilateral has vertices of all kinds (see Fig. 14 (right)). Let 14(1?) be the set of 



black vertices. A discrete S-isothermic surface is a map 

/b : Vb{V) ^ R\ 

with the following properties: 

(i) If wi, . . . , V2n S Vh{'D) are the neighbors of a ©-labeled vertex in cyclic order, 
then fb{vi),... , fb{v2n) lie on a circle in K'^ in the same cyclic order. This 
defines a map from the ©-labeled vertices to the set of circles in M'^. 
(a) If ui , . . . , V2n S VbCD) are the neighbors of an ©-labeled vertex, then fb{vi), . . . 
fb{v2n) lie on a sphere in K^. This defines a map from the ©-labeled vertices 
to the set of spheres in M.^. 
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(Hi) If Vc and Vg are the ©-labeled and the ©-labeled vertex of a quadrilateral of 
V, then the circle corresponding to Vc intersects the sphere corresponding to 
Vg orthogonally. 

Discrete S-isothermic surfaces are therefore composed of tangent spheres and 
tangent circles, with the spheres and circles intersecting orthogonally. The class of 
discrete S-isothermic surfaces is obviously invariant under Mobius transformations. 

Given a discrete S-isothermic surface, one can add the centers of the spheres 
and circles to it giving a map V{'D) R^. The discrete isothermic surface obtained 
is called the central extension of the discrete S-isothermic surface. All its faces are 
orthogonal kites. 

An important fact of the theory of isothermic surfaces (smooth and discrete) 
is the existence of a dual isothermic surface [6J. Let / : D D ^ M'' be an 
isothermic immersion. Then the formulas 

r* fx r* fy 

define an isothermic immersion f* : R"^ D D ^ which is called the dual 
isothermic surface. Indeed, one can easily check that the form df* is closed and 



/* satisfies (4.2 1. Exactly the same formulas can be applied in the discrete case. 



Proposition 4.4. Suppose f : T) ^ R is a discrete isothermic surface and suppose 



the edges have been consistently labeled '+ ' and as in Fig. 11 Then the dual 
discrete isothermic surface f*:'D^R^ is defined by the formula 

Af* = ± 

where Af denotes the difference of values at neighboring vertices and the sign is 
chosen according to the edge label. 

The closedness of the corresponding discrete form is elementary checked for 
one kite. 

Proposition 4.5. The dual of the central extension of a discrete S-isothermic surface 
is the central extension of another discrete S-isothermic surface. 

If we disregard the centers, we obtain the definition of the dual discrete S- 
isothermic surface. The dual discrete S-isothermic surface can be defined also 
without referring to the central extension 



5. Discrete Minimal Surfaces and Circle Patterns: Geometry from 
Combinatorics 

In this section (following 0]) we define discrete minimal S-isothermic surfaces 
(or discrete minimal surfaces for short) and present the most important facts from 
their theory. The main idea of the approach of [4j is the following. Minimal surfaces 
are characterized among (smooth) isothermic surfaces by the property that they 
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are dual to their Gauss map. The duality transformation and this characterization 
of minimal surfaces carries over to the discrete domain. The role of the Gauss map 
is played by discrete S-isothermic surfaces whose spheres all intersect one fixed 
sphere orthogonally. 

5.1. Koebe polyhedra 

A circle packing in 5^ is a configuration of disjoint discs which may touch but 
not intersect. The construction of discrete S-isothermic "round spheres" is based 
on their relation to circle packings in S^. The following theorem is central in this 
theory. 

Theorem 5.1. For every polytopa^cellular decomposition of the sphere, there exists 
a pattern of circles in the sphere with the following properties. There is a circle 
corresponding to each face and to each vertex. The vertex circles form a packing 
with two circles touching if and only if the corresponding vertices are adjacent. 
Likewise, the face circles form a packing with circles touching if and only if the 
corresponding faces are adjacent. For each edge, there is a pair of touching vertex 
circles and a pair of touching face circles. These pairs touch in the same point, 
intersecting each other orthogonally. 

This circle pattern is unique up to Mobius transformations. 

The first published statement and proof of this theorem is contained in |13j . 
For generalizations, see [33], [3T], and [9 , the latter also for a variational proof. 

Theorem |5.1| is a generalization of the following remarkable statement about 
circle packings due to Koebe ^5] . 

Theorem 5.2. (Koebe) For every triangulation of the sphere there is a packing of 
circles in the sphere such that circles correspond to vertices, and two circles touch 
if and only if the corresponding vertices are adjacent. This circle packing is unique 
up to Mobius transformations of the sphere. 

Observe that, for a triangulation, one automatically obtains not one but 



two orthogonally intersecting circle packings, as shown in Fig. 15 (right). Indeed, 
the circles passing through the points of contact of three mutually touching circles 
intersect these orthogonally, and thus Theorem |5.2| is a special case of Theorem |5.1| 

Consider a circle pattern of Theorem |5.1| Associating white vertices to circles 
and black vertices to their intersection points, one obtains a quad-graph. Actu- 
ally we have an S-quad-graph: Since the circle pattern is comprised by two circle 
packings intersecting orthogonally we have two kinds of white vertices, @ and ©, 
corresponding to the circles of the two packings. 

Now let us construct the spheres intersecting S*^ orthogonally along the circles 
marked by @. If we then connect the centers of touching spheres, we obtain a 
convex polyhedron, all of whose edges are tangent to the sphere S^. Moreover, the 



'^We call a cellular decomposition of a surface polytopal, if the closed cells are closed discs, and 
two closed cells intersect in one closed cell if at all. 
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Figure 15. Left: An orthogonal circle pattern corresponding to 
a cellular decomposition. Middle: A circle packing corresponding 
to a triangulation. Right: The orthogonal circles. 



circles marked with © are inscribed in the faces of the polyhedron. Thus we have 
a discrete S-isothermic surface. 



We arrive at the following theorem, which is equivalent to Theorem 5.1 (see 
also HO]). 

Theorem 5.3. Every polytopal cell decomposition of the sphere can be realized by 
a polyhedron with edges tangent to the sphere. This realization is unique up to 
projective transformations which fix the sphere. 

These polyhedra are called the Koebe polyhedra. We interpret the corre- 
sponding discrete S-isothermic surfaces as conformal discretizations of the "round 
sphere" . 

5.2. Definition of discrete minimal sm-faces 

Let / : P — > M'^ be a discrete S-isothermic surface. Suppose a; e I? is a white 
(s) vertex of the quad-graph I?, i.e. f{x) is the center of a sphere. Consider all 
quadrilaterals of T) incident to x and denote by t/i . . . y2n their black vertices and 
by xi . . . X2n their white © vertices. We will call the vertices xi . . . X2n neghbor- 
ing X in V. (Generically, n = 2.) Then f{yj) are the points of contact with the 
neighboring spheres and simultaneously points of intersection with the orthogonal 
circles centered at f{xj) (see Fig. 16 1. 



Consider only the white vertices of the quad-graph. Observe that each © 
vertex of a discrete S-isothermic surface and all its @ neighbors are coplanar. 
Indeed the plane of the circle centered at /(©) contains all its /(®) neighbors. 
The same condition imposed at the @ vertices leads to a special class of surfaces. 

Definition 5.4. A discrete S-isothermic surface / : I? ^ M'^ is called discrete mini- 
mal if for each sphere of /(2?) there exists a plane which contains the center /(©) 
of the sphere as well as the centers /(©) of all neighboring circles orthogonal to the 
sphere, i.e. if each white vertex of /(2?) is coplanar to all its white neighbors. These 
planes should be considered as tangent planes to the discrete minimal surface. 

Theorem 5.5. An S-isothermic discrete surface f is a discrete minimal surface, if 
and only if the dual S-isothermic surface f* corresponds to a Koebe polyhedron. 
In that case the dual surface N := f* : Vt^(T>) — > M.^ at white vertices V^(T>) can 
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Figure 16. To the definition of discrete minimal surfaces: the 
tangent plane through the center f{x) of a sphere and the cen- 
ters f{xj) of the neighboring circles. The circles and the sphere 
intersect orthogonally at the black points fiuj)- 



be treated as the Gauss map N of the discrete minimal surface: At (s) vertices N is 
orthogonal to the tangent planes, and at ©vertices N is orthogonal to the planes 
of the circles centered at /(©). 

Proof. That the S-isothermic dual of a Koebe polyhedron is a discrete minimal 
surface is fairly obvious. On the other hand, let / : I? ^ E'^ be a discrete minimal 
surface and let x G V and yi ■ ■ ■ y2n G I? be as above. Let /* : I? — > K'^ be the 
dual S-isothermic surface. We need to show that all circles of /* lie in one and the 
same sphere S and that all the spheres of /* intersect S orthogonally. Since the 
quadrilaterals of a discrete S-isothermic surface are kites the minimality condition 



of Definition 5.4 can be reformulated as follows: There is a plane through f{x) 
such that the points {fiyj) \ j even} and the points {f{yj) \ j odd} lie in planes 



which are parallel to it at the same distance on opposite sides (see Fig. 16 1. It 
follows immediately that the points /*(yi) ■ • . /*(2/2n) lie on a circle in a sphere 
Sx around f*{x). The plane of Cx is orthogonal to the normal N to the tangent 
plane at f{x). Let S be the sphere which intersects Sx orthogonally in Cx- The 
orthogonal circles through f*{yi) . . . f*{y2n) also lie in S. Hence, all spheres of /* 
intersect S orthogonally and all circles of /* lie in S. □ 



Theorem 5.5 is a complete analogue of Christoffel's characterization [15 of 



continuous minimal surfaces. 

Theorem 5.6. (Christoffel) Minimal surfaces are isothermic. An isothermic im- 
mersion is a minimal surface, if and and only if the dual immersion is contained 
in a sphere. In that case the dual immersion is in fact the Gauss map of the 
minimal surface. 
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Thus a discrete minimal surface is a discrete S-isothermic surface which is 
dual to a Koebe polyhedron; the latter is its Gauss map and is a discrete analogue 
of a conformal parametrization of the sphere. 

The simplest infinite orthogonal circle pattern in the plane consists of circles 
with equal radius r and centers on a square grid with spacing ^^/2r. One can 
project it stereographically to the sphere, construct orthogonal spheres through 
half of the circles and dualize to obtain a discrete version of Enneper's surface. See 
Fig. [Tlfri^/ifj. Only the circles are shown. 

5.3. Construction of discrete minimal sm-faces 

A general method to construct discrete minimal surfaces is schematically shown 
in the following diagram. 

continuous minimal surface 
image of curvature lines under Gauss-map 
cell decomposition of (a branched cover of) the sphere 
orthogonal circle pattern 
Koebe polyhedron 

discrete minimal surface 

As usual in the theory on minimal surfaces [25j . one starts constructing such a 
surface with a rough idea of how it should look. To use our method, one should 
understand its Gauss map and the combinatorics of the curvature-line pattern. 
The image of the curvature-line pattern under the Gauss map provides us with 
a cell decomposition of (a part of) S'^ or a covering. From these data, applying 
Theorem |5.1| we obtain a Koebe polyhedron with the prescribed combinatorics. 
Finally, the dualization step yields the desired discrete minimal surface. For the 



discrete Schwarz P-surface the construction method is demonstrated in Fig. 17 
and Fig. [TS] 

Let us emphasize that our data, aside from possible boundary conditions, are 
purely combinatorial — the combinatorics of the curvature line pattern. All faces are 
quadrilaterals and typical vertices have four edges. There may exist distinguished 
vertices (corresponding to the ends or umbilic points of a minimal surface) with a 
different number of edges. 

The most nontrivial step in the above construction is the third one listed 
in the diagram. It is based on the generalized Koebe Theorem |5.1| It implies 
the existence and uniqueness for the discrete minimal S-isothermic surface under 
consideration, but not only this. A constructive proof of the generalized Koebe 
theorem suggested in [9J is based on a variational principle and also provides a 
method for the numerical construction of circle patterns. An alternative algorithm 
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Gauss image of the curvature lines circle pattern 




Koebe polyhedron discrete minimal surface 



Figure 17. Construction of the discrete Schwarz P-surface. 

by Thurston was implemented in Stephenson's program circlepack (see |38j for 
an exhaustive presentation of the theory of circle packings and their numerics) . The 
first step is to transfer the problem from the sphere to the plane by a stereographic 
projection. Then the radii of the circles are calculated. If the radii are known, it 
is easy to reconstruct the circle pattern. The above-mentioned variational method 
is based on the observation that the equations for the radii are the equations for 
a critical point of a convex function of the radii. The variational method involves 
minimizing this function to solve the equations. 

Let us describe the variational method of [9 for construction of (orthogonal) 
circle patterns in the plane and demonstrate how it can be applied to construct 
the discrete Schwarz-P surface. 

Instead of the radii r of the circles, we use the logarithmic radii 

p — log r. 

For each circle j, we need to find a pj such that the corresponding radii solve the 
circle pattern problem. This leads to the following equations, one for each circle. 
The equation for circle j is 



2 ^ arctan e^"' = <i>j , 

neighbors k 



(5.1) 



Surfaces from Circles 



27 



where the sum is taken over all neighboring circles k. For each circle j, <i>j is the 
nominal angle covered by the neighboring circles. It is normally 27r for interior 
circles, but it differs for circles on the boundary or for circles where the pattern 
branches. 

Theorem 5.7. The critical points of the functional 
Sip) = J2 {l^U2iieP''-P')+linU2iieP'''"') - ^(p, + Pfc)) + ^ 

(j,k) j 

correspond to orthogonal circle patterns in the plane with cone angles at the 
centers of circles (^j = 27r for internal circles). Here, the first sum is taken over 
all pairs (j, k) of neighboring circles, the second sum is taken over all circles j, and 
the dilogarithm function L\2{z) is defined by Li2(z) = — JJ^log(l — C)dC/C- The 
functional is scale-invariant and restricted to the subspace pj — it is convex. 



Proof. The formula for the functional follows from (5.1 ) and 



d 1 1 + ie^ 

— ImLi2(ie^) = — log = arctane"^. 

dx 2i 1 — ie^ 

The second derivative of the functional is the quadratic form 

= E —u^ — ^ (^/^fe - ^^^)'' 

(^)Cosh(pfc-pj) 

where the sum is taken over pairs of neighboring circles, which implies the con- 
vexity. □ 

Now the idea is to minimize S{p) restricted to pj — 0. The convexity of 
the functional implies the existence and uniqueness of solutions of the classical 
boundary valued problems: Dirichlet (with prescribed pj on the boundary) and 
Neumann (with prescribed <i>j on the boundary). 

The Schwarz-P surface is a triply periodic minimal surface. It is the symmetric 
case in a 2-parameter family of minimal surfaces with 3 different hole sizes (only 
the ratios of the hole sizes matter), see [T7]. The Gauss map is a double cover of 
the sphere with 8 branch points. The image of the curvature line pattern under the 
Gauss map is shown schematically in Fig. [Tt] (top left), thin lines. It is a refined 
cube. More generally, one may consider three different numbers m, n, and k of 
slices in the three directions. The 8 corners of the cube correspond to the branch 
points of the Gauss map. Hence, not 3 but 6 edges are incident with each corner 
vertex. The corner vertices are assigned the label ©. We assume that the numbers 
TO, n, and k are even, so that the vertices of the quad graph may be labelled '©', 

and '•' consistently (see Section|4|. 

We will take advantage of the symmetry of the surface and construct only a 
piece of the corresponding circle pattern. Indeed the combinatorics of the quad- 
graph has the symmetry of a rectangular papallelepipedon. We construct an or- 
thogonal circle pattern with the symmetry of the rectangular papallelepipedon 
eliminating the Mobius ambiguity of Theorem |5.1[ Consider one fourth of the 
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sphere bounded by two orthogonal great circles connecting the north and the 
south poles of the sphere. There are two distinguished vertices (corners of the 
cube) in this piece. Mapping the north pole to the origin and the south pole to 
infinity by stereographic projection we obtain a Neumann boundary value problem 
for orthogonal circle patterns in the plane. The symmetry great circles become two 
orthogonal straight lines. The solution of this problem is shown in Fig. [19] The 
Neumann boundary data are $ = 7r/2 for the lower left and upper right boundary 
circles and $ = tt for all other boundary circles (along the symmetry lines). 

Now map this circle pattern to the sphere by the same stereographic projec- 
tion. One half of the spherical circle pattern obtained (above the equator) is shown 



in Fig. 17 (top right). This is one eighth of the complete spherical pattern. Now lift 



the circle pattern to the branched cover, construct the Koebe polyhedron and du- 



alize it to obtain the Schwarz-P surface; see Fig. 17 (bottom row). A translational 
fundamental piece of the surface is shown in Fig. |18[ 
We summarize these results in a theorem. 

Theorem 5.8. Given three even positive integers m, n, k, there exists a correspond- 
ing unique (unsymmetric) S-isothermic Schwarz-P surface. 

Surfaces with the same ratios m : n : k are diS^erent discretizations of the 
same continuous Schwarz-P surface. The cases with m = n = k correspond to the 
symmetric Schwarz-P surface. 

Further examples of discrete minimal surfaces can be found in |4] and [14] . 

Orthogonal circle patterns on the sphere are treated as the Gauss map of the 
discrete minimal surface and are central in this theory. Although circle patterns 
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Figure 19. A piece of the circle pattern for a Schwarz-P surface 
after stereographic projection to the plane. 



on the plane and on the sphere differ just by the stereographic projection some 
geometric properties of the Gauss map can get lost when represented in the plane. 
Moreover to produce branched circle patterns in the sphere it is important to be 
able to work with circle patterns directly on the sphere. A variational method which 
works directly on the sphere was suggested in |35L ^ . This variational principle 
for spherical circle patterns is analogous to the variational principles for Euclidean 
and hyperbolic patterns presented in [S]. Unlike the Euclidean and hyperbolic 
cases the spherical functional is not convex, which makes difficult to use it in the 
theory. However the spherical functional has proved to be amazingly powerful for 
practical computation of spherical circle patterns (see [3 for detail). In particular, 
it can be used to produce branched circle patterns in the sphere. 

Numerous examples of discrete minimal surfaces are constructed with the 
help of the spherical functional in the contribution of Biicking |14j to this volume. 
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6. Discrete Conformal Surfaces and Circle Patterns 



Conformal immersions 



f -.R^ D D 




fx-^fy ■ 



Any surface can be conformally parametrized, and a conformal parametrization 
induces a complex structure on D in which z — x + iy is a, local complex coor- 
dinate. The development of a theory of discrete conformal meshes and discrete 
Riemann surfaces is one of the popular topics in discrete differential geometry. 
Due to their almost square quadrilateral faces and general applicability, discrete 
conformal parametrizations are important in computer graphics, in particular for 
texture mapping. Recent progress in this area could be a topic of another paper; 
it lies beyond the scope of this survey. In this section we mention shortly only the 
methods based on circle patterns. 

Conformal mappings can be characterized as mapping infinitesimal circles to 
infinitesimal circles. The idea to replace infinitesimal circles with finite circles is 
quite natural and leads to circle packings and circle patterns (see also Section |4| . 
The corresponding theory for conformal mappings to the plane is well developed 
(see the recent book of Stephenson [3H]). The discrete conformal mappings are 
constructed using a version of Koebe's theorem or a variational principle which 
imply the corresponding existence and uniqueness statements as well as a numeri- 
cal method for construction (see Section [s]). It is proven that a conformal mapping 
can be approximated by a sequence of increasingly fine, regular circle packings or 
patterns [Ml^ . 

Several attempts have been made to generalize this theory for discrete con- 
formal parametrizations of surfaces. 

The simplest natural idea is to ignore the geometry and take only the com- 
binatorial data of the simplicial surface |38]. Due to Koebe's theorem there exists 
essentially unique circle packing representing this combinatorics. One treats it as 
a discrete conformal mapping of the surface. This method has been successfully 
applied by Hurdal et al. for fiat mapping of the human cerebellum. However 
a serious disadvantage is that the results depend only on the combinatorics and 
not on the geometry of the original mesh. 

An extension of Stephenson's circle packing scheme which takes the geometry 
into account is due to Bowers and Hurdal |12j . They treat circle patterns with non- 
intersecting circles corresponding to vertices of the original mesh. The geometric 
data in this case are the so called inversive distances of pairs of neighboring circles, 
which can be treated as imaginary intersection angles of circles. The idea is to get 
a discrete conformal mapping as a circle pattern in the plane with the inversive 
distances coinciding with the inversive distances of small spheres in space. The 
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latter are centered at the vertices of the original meshes. The disadvantage of this 
method is that there are almost no theoretical results regarding the existence and 
uniqueness of inversive distance patterns. 

Kharevich, Springborn and Schroder suggested another way ;27J to handle the 
geometric data. They consider the circumcircles of the faces of a simplicial surface 
and take their intersection angles 9. The circumcircles are taken intrinsically, i.e. 
they are round circles with respect to the surface metric. The latter is a flat 
metric with conical singularities at vertices of the mesh. The intersection angles 
9 are the geometric data to be respected, i.e. ideally for a discrete conformal 
mapping one wishes a circle pattern in the plane with the same intersection angles 
9. An advantage of this method is that similarly to the circle packing method it is 
based on a solid theoretical background — the variational principle for patterns of 
intersecting circles [9 . However to get a circle pattern flat one has to change the 
intersection angles 9 9 and it seems that there is no geometric way to do this 
(in [27 the angles 9 of the circle pattern in the plane are defined as minimizing 
the sum of differences squared {9 — 9)^). A solution to this problem could possibly 
be achieved by a method based on Delaunay triangulations of circle patterns with 
disjoint circles. The corresponding variational principle has been found recently in 
[36]. 

It seems that discrete conformal surface parametrizations are at the beginning 
of a promising development. Although now only some basic ideas about discrete 
conformal surface parametrizations have been clarifled and no approximation re- 
sults are known, there is good chance for a fundamental theory with practical 
applications in this field. 
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